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The distinct distribution of local magnetic fields due to superconducting vortices can be detected 
with nuclear magnetic resonance (NMR) and used to investigate vortices and related physical prop- 
erties of extreme type II superconductivity. This review summarizes work on high temperature 
superconductors (HTS) including cuprates and pnictide materials. Recent experimental results are 
presented which reveal the nature of vortex matter and novel electronic states. For example, the 
NMR spectrum has been found to provide a sharp indication of the vortex melting transition. In 
the vortex solid a frequency dependent spin-lattice relaxation has been reported in cuprates, in- 
cluding YBa2Cu307_a;, Bi2SrCa2Cu208+i , and Tl2Ba2Cu06+5 . These results have initiated a new 
spectroscopy via Doppler shifted nodal quasiparticles for the investigation of vortices. At very high 
magnetic fields this approach is a promising method for the study of vortex core excitations. These 
measurements have been used to quantify an induced spin density wave near the vortex cores in 
Bi2SrCa2Cu208+i . Although the cuprates have a different superconducting order parameter than 
the iron arsenide superconductors there are, nonetheless, some striking similarities between them 
regarding vortex dynamics and frequency dependent relaxation. 



PACS numbers: 74.25.Uv, 74.72.-h, 75.30.Fv 
INTRODUCTION 



Among the many tools used by experimenters to ex- 
plore the properties of superconductors NMR has played 
an important role. In this brief review we focus on appli- 
cations of nuclear magnetic resonance (NMR) to study 
the mixed state in type II superconductivity. The in- 
terest in this topic and the corresponding explosion in 
the number of publications where NMR has been of cen- 
tral importance, came immediately after the discovery of 
high temperature superconductivity^] with the earliest 
NMR and nuclear quadrupole resonance (NQR) papers 
reporting on (Lai_a:Sra;)2Cu04,[2] YBa2Cu307_a;,t3 and 
Tl2Ba2Ca2Cu3O10-1-5.ll! NMR is a microscopic probe of 
the electronic state and has the advantage of being spec- 
troscopically selective for specific crystallographic posi- 
tions in the superconducting compound. The spectrum 
is the distribution of local magnetic fields for which there 
are several contributions from superconductivity. There 
are fields from diamagnetic currents circulating at the 
surface of the sample, especially important at low ap- 
plied magnetic fields, superposed on the inhomogeneous 
distributions of magnetic field from vortex supercurrents 
from the sample interior. An ideal distribution, calcu- 
lated from GL theory using algorithms from Brandt, [3 [6] 
is shown in Fig. 1 and 2, called the Redfield pattern. Ad- 
ditionally, the spin-lattice relaxation rate, T^^, gives a 
measure of the electronic excitations which provide im- 
portant signatures of the superconducting state. These 
aspects have been described in the early literature, no- 
tably in the review by MacLaughlin[7] that emphasizes 
the importance of NMR in providing early support for 
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FIG. 1. The vortex distribution of magnetic fields equiva- 
lent to the NMR spectrum often called the Redfield pattern. 
The singularities in the distribution are labeled: the minimum 
field, m; the saddle point field, s; and the maximum field, M, 
at the vortex core. The inset shows the corresponding spa- 
tial distribution. This figure was taken from Mitrovic et al. [S] 
calculated from Brandt's algorithm 5 for // = 37 T. 

the BCS theory. 

After the discovery of high temperature superconduc- 
tivity (HTS), a large number of excellent reviews of NMR 
have been written; among these are: Pennington and 
Slichter,[3] Asayama et al., [TU] Berthier et al., [TT] Riga- 
monti et al, [12^ Walstedt,[T3] and Curro. [M] This work 
has been followed by the discovery of new compounds, 
or new work on previously known materials, that have 
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challenged our understanding of the symmetry of the or- 
der parameter in the superconducting state, notably for 
UPta, Sr2Ru04, CeCoIns, MgB2, cuprates, and the iron 
pnictides and their related compounds. In this review we 
will consider aspects of vortex structures in a selection of 
these materials that have been recently explored using 
NMR. Two important advances in the past decade in ap- 
plications of NMR to the understanding of vortex struc- 
tures in HTS are that high quality single crystals have 
become available and measurements at very high mag- 
netic fields have become routine. Most of the focus in 
this review will be directed toward these developments. 
Contributions to our survey are largely taken from the 
condensed matter NMR group at the National High Mag- 
netic Field Laboratory (NHMFL) in Tallahassee, Florida, 
where many of these developments were made, and their 
collaborations with the NMR group at Northwestern Uni- 
versity. 

VORTEX SPECTRUM 

It was shown by Abrikosov|15| for type II supercon- 
ductors, K = A/^ >> 1, that quantized vortices pene- 
trate a superconductor with a penetration depth. A, and 
coherence length ^. The Ginzburg-Landau parameter, 
K ^ 100, satisfies this condition for HTS. From the Lon- 
don equations, |16] it can be shown that the magnetic 
field distribution takes the form shown in Fig. 1 and 2 
where the quantized fiux bundles, or fiuxons, are assumed 
to be rectilinear. In this London model for the vortex 
structure the vortex core region is excluded since it is 
assumed that the coherence length is sufficiently small 
and the external magnetic field, H , is substantially lower 
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FIG. 2. The vortex spectrum, supercurrent momentum, pa, 
supercurrent density, js, and order parameter, OP, calculated 
by Mitrovic'6 for YBa2Cu307-.a; using Brandt's method^' 
at H = 37T.' 



than the upper critical field Hc2 ■ From Ginzburg-Landau 
theory O [15] it was shown that the vortices form a lattice 
and within this structure there is a spatially inhomoge- 
neous distribution of magnetic fields, supercurrent mo- 
menta, supercurrent density, and the order parameter. 
Using algorithms developed by Brandt [S] these distribu- 
tions can be calculated for any specific values of A and ^, 
as is shown in Fig. 2 taken from Mitrovic's calculation [6] 
for YBa2Cu307_a;. The majority of type II supercon- 
ductors are highly anisotropic with an anisotropy axis 
parallel to the crystal c-axis. The flux distribution for 
H\\c can be expected to form a two-dimensional lattice 
with a high degree of symmetry, such as square or hexag- 
onal, or simple distortions of these. Such details will alter 
the Redfield pattern, but not sufficiently to be of concern 
here. In fact the properties of the inhomogeneous mixed 
state are less dependent on this symmetry approaching 
the vortex core and at high magnetic fields. 

The classic Redfield pattern from the GL-theory 
is not generally found to be in very good quanti- 
tative agreement with NMR experiments. Usually 
there are broadening mechanisms in addition to the 
local fields from supercurrents which are convolved 
with the Redfield pattern, as for example, in the 
works on YBaaCuaOr-a;.!!! EZ] Tl2Sr2BaCu206.8,lIS] 
Tl2Ba2Cu06+5,ll9j and Bi2SrCa2Cu2 08+5.^20] An ex- 
ample where good agreement was obtained is for the su- 
perconducting state of CeCoIns shown in Fig. 3. Another 
case is the lineshape simulation of the powder pattern 
of the strongly anisotropic superconductor, MgB9.pT] In 
this instance the frequency shift from the magnetic mo- 
ment of the screening currents was explicitly included in 
the GL-calculation along with the corresponding mag- 
netic shift distribution from the vortex supercurrents. 
An excellent representation of the full, quadrupolar-split, 
spectrum was obtained. In these instances the penetra- 
tion depth, and with less accuracy the coherence length, 
can be inferred as fitting parameters. However, without 
a consistent representation of the lineshape from the the- 
ory a determination of the penetration depth from either 
the second moment of the spectrum or the full-width-at- 
half- maximum is not very meaningful. [22] 

VORTEX PHASES 

For NMR to probe the vortex solid structure, vortices 
must be stationary on the time scale of the NMR exper- 
iment. As thermal fluctuations increase at higher tem- 
peratures, vortices enter a liquid-like phase averaging out 
the inhomogeneous magnetic fields and motionally nar- 
rowing the NMR spectrum. For low transition temper- 
ature superconductors, the thermal fluctuations are less 
energetic resulting in a robust solid vortex phase, stable 
up to Tc- However, in the case of HTS, with its high tran- 
sition temperatures and high anisotropy, the result is a 
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FIG. 3. The vortex distribution of magnetic fields in CeCoIns 
taken from Koutroulakis et al.[23 measured from their ^^^In 
NMR spectrum and compared with their calculated spectrum 
from Brandt's algorithm|5j at = 4 T. 
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FIG. 4. The vortex phase diagram for Bi2SrCa2Cu208+i 
taken from Chen et al.^7\ According to the theory of Glaz- 
man and Koshelev, [31] there is a limiting value of the melting 
transition at high magnetic field (solid red line) which, theo- 
retically, should be an ideal two-dimensional transition. From 
experiment it was found that T2D = 12 K for an overdoped 
crystal with Tc = 75 K and corresponds to a penetration 
depth of Xab = 220 nm. The dotted line is a fit to the theory. 



stable vortex liquid that solidifies only at sufficiently low 
temperatures. In some cases, such as highly anisotropic 
Bi2SrCa2Cu208+5, the liquid phase can extend over a 
significant region of the magnetic field-temperature phase 
diagram. Fig. [4j Vortex melting is in fact the only true 
thermodynamic transition for a clean type II supercon- 
ductor in a magnetic field. [21] The first order signature 
of this transition has been extensively studied for both 
YBa2Cu307-:E[25; and Bi2SrCa2Cu2 08+5[3i H?]. It is 
essential that NMR experiments intended to explore spa- 
tial configurations of vortices be performed well into the 
vortex solid domain. 

Formation of a vortex solid has been identified from 
NMR measuring either spin-spin relaxation, T2^, or 
the NMR linewidth as a function of temperature. For 
measurements on YBa2Cu307_x an abrupt change in 
^^O T^^ was observed [25] and related to a change in 
vortex dynamics. This corresponded well to the irre- 
versibility line at least at low magnetic fields. Sim- 
ilar behavior, even more pronounced, has been ob- 
served in ^^As Tj"^ experiments on the electron doped 
pnictide, Ba(Feo.93Coo.o7)2As2.[29j Additionally, changes 
in the linewidth have been helpful to identify vor- 
tex freezing as reported for YBa2Cu307_a; [281 IHO] and 
Ba(Feo.93Coo.o7)2As2.Bg 

In the case of ^^O NMR on Bi2SrCa2Cu208+5 there 
is an abrupt signature of the formation of the solid 
vortex phase at a well defined temperature where the 
magnetic field distribution increases abruptly on cool- 
ing, determined from the second moment of the NMR 
spectrum. |27j In Fig. 4 we show the measurement of 



the phase diagram which results from plotting this sig- 
nature, compared to two-dimensional melting theorv[5T| 
with which there is good agreement at high magnetic 
field. The theory shows that the vortex-vortex inter- 
actions are largely determined by the electromagnetic 
coupling between supercurrents; the contributions from 
Josephson tunneling between cuprate layers are relatively 
suppressed in the limit of high field leading to quasi two- 
dimensional behavior. The phase boundary was calcu- 
lated to approach a vertical asymptote for the limit of 
an ideal two-dimensional melting transition which corre- 
sponds well to experimental observations. j27j 



SPATIALLY RESOLVED NMR 

As it became more evident that the high Tc super- 
conducting cuprates have an order parameter with d- 
wave symmetry, it was recognized that the existence 
of nodes in the gap at the Fermi surface would affect 
the thermodynamics at low temperatures owing to vor- 
tex supercurrents. [32] The spatially averaged density of 
states near the Fermi surface was calculated for d-wave 
pairing symmetry [33] including Doppler contributions to 
the quasiparticle excitation spectrum from the vortex 
supercurrents. This term in the energy has the form, 
Se ^ vf Ps, where Ps oc l/r, and is the Fermi ve- 
locity, Ps is the supercurrent momentum, and r is the 
distance from the vortex core. It was found that the spa- 
tially averaged density of states has a unique ^/H depen- 
dence, which became known as the Volovik effect. Bulk 
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FIG. 5. The magnetic field dependence of the "^^Cu NMR re- 
laxation rate is shown for Tl2Sr2BaCu206.8 taken from Zheng 
et al. [18] .The linear behavior is expected for the average rate 
according to the Volovik effect for a d-wave superconductor. 
There is no evidence here for a Zeeman contribution to the 
quasiparticle excitations which is expected[8] to vary with 
magnetic field according to oc H^. 

experimental probes of the density of states including 
specific heat [Ml [55] and thermal transport [SS] [37] were 
in agreement with this picture. 

Similarly, the average of T^^ should be proportional to 
the average of the joint density of states that is expected 
to depend on magnetic field, oc H. This behavior has 
been reported in Tl2Sr2BaCu206.8 by Zheng et aZ.[T5] 
and is displayed in Fig. [5] Further discussion of the 
Volovik effect is deferred to appendix A. 

However, the behavior of the vortex core and extended 
quasiparticle states on a local microscopic scale were still 
not well-understood and the need for a local probe of 
the electronic excitations was evident. Scanning tunnel- 
ing microscopy and scanning tunneling spectroscopy are 
sensitive to the local density of states; but, these probes 
are limited to the surface of the sample in contrast to 
NMR. The sample volume to which NMR responds is 
within a London penetration depth. A, of the surface, 
generally a few hundred nm in size. Superconductivity is 
suppressed very close to the surface of the sample, on the 
much smaller length scale given by the superconducting 
coherence length, ^, typically a few nm for extreme type 
II materials. This leaves a region of the superconductor 
to which NMR is sensitive, near the surface but suffi- 
ciently far away that it is unperturbed by surface effects. 

It was predicted theoretically that NMR relaxation 
would be sensitive to spatially dependent quasiparticle 
states, [M] since the spin-lattice relaxation rate, Tj"^, de- 
pends on the density of states, N{e). As pointed out 
above, in the mixed state of the superconductor, N(e) 



must vary in space in response to Doppler contributions 
from the vortex supercurrents. With a well-defined vor- 
tex lattice spectrum, the local magnetic fields correspond 
to points in real space, see the inset to Fig. 1. So it 
was suggested that examining the relaxation at different 
internal magnetic fields, resolved within the NMR spec- 
trum at different freqiencies, would allow a determination 
of quasiparticle excitations to be mapped throughout real 
space in what we call spatially resolved NMR.|5] In the 
following we will discuss various examples of spatially re- 
solved NMR experiments on HTS. 

Spin-Lattice Relaxation: Experiment 

The NMR relaxation rate in a normal metal is calcu- 
lated in terms of the thermal average of the joint density 
of states: 

l/Ti cx J iV(6,)A^(e/)/(e,)[l - /(e/)] de (1) 

where i and / indices label the initial and final elec- 
tron states which are necessarily spin-up and spin-down 
in order to cause a nuclear spin flip and conserve to- 
tal spin angular momentum, and /(e) is the Fermi-Dirac 
distribution function. The change in energy of the nu- 
clear spin is essentially zero so = e/, but the excita- 
tions can come from nodes at different positions, k^j, on 
the Fermi surface. Arguments extending this framework 
to the superconducting state involve coherence factors [7] 
important in elemental BCS superconductors, but which 
have been irrelevant to the class of HTS superconductors, 
largely owing to strong anisotropy. Furthermore, for ex- 
otic or unconventional pairing systems there is an inter- 
play between order parameter structure, impurity scat- 
tering, and high magnetic fields, together with the shift 
of the local electronic energy spectrum from the Doppler 
effect, which requires sophisticated analysis. One limit 
that is arguably tractable is that of low temperatures for 
nodal superconductors where the excitations are Dirac- 
like. In this case the thermal contributions to the integral 
in Eq. 1 are less important and limiting forms of A''(e) 
at low energy can be used. For cuprates this is simply, 
N{e) oc e. 

It is a common error in interpretation of T^^ in the 
superconducting state to look for power law temperature 
dependences below Tc, say in the region ~0.5 < T/Tc < 
1.0 where there is a sudden drop in the relaxation rate 
compared to the normal state of the metal. This type 
of analysis corresponds to an inappropriate expansion of 
Eq. 1 in a region in which the structure of the order pa- 
rameter has not approached its low temperature limiting 
form. Impurity effects on N{e) are difficult to quantify 
and, in principle, cannot be excluded even for nominally 
clean single crystals of superconductors. However, the 
low temperature limit is potentially accessible in clean 
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samples at sufficiently high magnetic fields such that the 
field energy is larger than the impurity bandwidth. 

Experimental evidence for spatially dependent relax- 
ation came from ^^O NMR in isotopically enriched, op- 
timally doped, ahgned powders of YBa^Cu^OT-^.p 
In the first experiments by Curro et a/., [33] the central 
transition of the planar oxygen was isolated from the 
overlapping apical oxygen resonance by a technique of 
population-enhanced, double resonance. |40j The relax- 
ation rate was measured at different points in the NMR 
spectrum in a magnetic field oi H = 8.3 T at several 
temperatures below the vortex melting transition tem- 
perature. The relaxation was found to be inhomoge- 
neously distributed and increasing with Hint, i-e. with 
frequency within the NMR spectrum for the planar ^^O 
site. At temperatures below Tc the apical T^^ has a T- 
linear temperature dependence, whereas the planar T^^ 
is cx T"^ upon entering the superconducting phase, but 
becomes T^^ oc T at low temperatures. These authors 
also measured and both rates were interpreted as 

evidence of vortex vibrations for T < 25 K and from 
quasiparticle scattering for T > 25 in the vortex solid 
state. However, theoretical investigations of the effect of 
vortex vibrations on T^^ by Wortis et oZ. [41] have sug- 
gested that this interpretation is unlikely. 

Mitrovic et al.^ [33] reported results on isotopi- 
cally enriched ^^O optimally-doped, aligned powders of 
YBa2Cu307_2:. Their work was performed up to very 
high magnetic fields, H < 42 T, Fig. [6ja). These au- 
thors came to a different conclusion than Curro et al. [39] 
regarding the low temperature mechanism for relaxation. 
Mitrovic et al. examined the -1/2 o -3/2 transition to 
separate the planar and apical oxygen resonances. The 
temperature and magnetic field dependence of T^^ re- 
laxation was analyzed in terms of scattering of quasipar- 
ticles between the gap nodes on the Fermi surface. The 
key features of their experiments are the dependence of 
T^^ on magnetic field and extension of measurements to 
very high fields to permit identification of the relaxation 
mechanism and sensitivity to vortex core excitations. At 
low temperatures the NMR spin-lattice relaxation rate 
due to quasiparticle scattering is given by, [5] 

T^;^ ^ {\e- Z + D,\\e + Z + Df\) (2) 

where e ~ ksT is the thermal contribution, Z — 
is the Zeeman energy, and Dij = {vp)ij-ps are the ini- 
tial and final Doppler shifts. By changing the external 
magnetic field, temperature, and position in the NMR 
spectrum, each of these terms changes magnitude. 

With increasing magnetic field, T^^ increases quadrat- 
ically with H at the saddle point of the spectrum where 
D is minimal. In this scenario Tj~^~ — e^|, and at 
low enough temperature becomes ~ \Z'^\, since the ther- 
mal term varies as T^. This behavior is demonstrated 
in Fig. [7] that compares the behaviors of YBa2Cu307_i: 
and Bi2SrCa2Cu208+5 showing an increase in T^^ cx 



H^. At locations where the Doppler term is more preva- 
lent and at low temperatures, T^^ is given by \ ±D'^ — Z^\ 
and the sign of D depends on the nodes for initial and 
final quasiparticle states. Observation of an increasing 
Tj~^ with magnetic field suggests [Sj that internode scat- 
tering dominates, ^ k/. 

Experiments by Kakuyanagi et al, 42J Fig. |6](b), re- 
port similar effects. The authors concentrate on the 
non-monotonic variations of T^^ across the spectrum. 
Separating the planar and apical resonances by exam- 
ining the -1/2 o -3/2 transition, they found different 
behavior at the different oxygen lattice sites. The pla- 
nar oxygen shows an increase in T^^ across the spec- 
trum which reaches a maximum and then diminishes at 
high frequency, approaching the vortex core. The api- 
cal oxygen shows a more slowly increasing T^^ , with a 
magnitude 5 times less than the planar site. The differ- 
ent behavior and magnitude of T^^ at these sites seems 
to indicate a different mechanism between the two loca- 
tions, specifically these authors conclude that the planar 
oxygen T^^ is dominated by quasiparticle interactions. 
They argue that if there were to be vortex vibrations they 
would affect the two sites in the same way. The dip in 
^ at high frequency, a characteristic inconsistent with 
Doppler shift, was interpreted as an effect of the vortex 
core, where superconductivity is suppressed, Fig. [6](b). 
Also, an increase on the low frequency side of the spec- 
trum could not be explained in terms of Doppler effects. 
There is a striking similarity in the two measurements of 
Tj"^ on YBa2Cu307_a: and YBa2Cu408 for comparable 
magnetic fields, 7J < 10 T. 

This body of experimental work paints a reasonably 
consistent picture. Certainly T^^ increases across 
the NMR spectrum due to vortices, and another 
contribution other than the Doppler shift affects 
T^^ near the vortex core, corresponding to the high 
frequency tail of the spectrum. The small upturn 
in Tj"^ on the low frequency part of the spectra is 
a puzzle, and appears to be ubiquitous to all re- 
ports to date: YBa2Cu307-:E,[i HI] YBa2Cu408,[31] 
Tl2Sr2BaCu206.8,[IH] Tl2Ba2Cu06+5,[IS] and for 
Bi2SrCa2Cu208+5.[2n] However, theoretical input helps 
clarify possible mechanisms for the inhomogeneous 
relaxation and further detailed experiments in the 
vortex core region of the spectrum will be essential to 
reveal the spin character of the vortex cores in HTS. 

Spin-Lattice Relaxation: Theory 

Early experimental indications of spatially inhomoge- 
neous relaxation in YBa2Cu307_:r quickly spurred sev- 
eral theoretical models. The two T^^ mechanisms are 
the effect of vortex vibrations and the spin flip scattering 
of quasiparticles. Wortis et al. [H] compared the relax- 
ation due to these two processes for ^'^Cu. Their cal- 
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FIG. 6. (a) The frequency dependent ^^O relaxation at various magnetic fields in YBa2Cu307_a; from Mitrovic et aZ. [8] ; 
and (b) for YBa2Cu408 at various temperatures from Kakuyanagi et a/.jj^ For YBa2Cu307_a;, panel (a), as the external 
magnetic field is increased the relaxation rate at each frequency position in the spectrum shows a constant increase attributed 
to the Zeeman effect. Additionally, T^^ increases with frequency across the spectrum consistent with a Doppler effect. The 
shaded region corresponds to the fraction of the measured ^'^O spectrum, shown in the background, that would correspond 
to nuclei in the vortex core for H = 37 T. For panel (b), when YBa2Cu408 is cooled through the superconducting transition, 
the relaxation shows no frequency dependence as would be expected for a motionally averaged vortex liquid. However, as the 
sample is cooled below the vortex freezing temperature, the relaxation becomes frequency dependent, characteristic of Doppler 
shifted relaxation with the exception of the highest and lowest frequency portions of the spectrum. The two experiments, (a) 
and (b), have clear similarities. 



culations resulted in spatially inhomogeneous relaxation 
which increased with internal magnetic field for both 
mechanisms, however the calculated T^^ due to vortex 
vibrations is much too slow to account for experimental 
results. 

There are two regimes examined for spin flip scatter- 
ing of quasiparticles; the weakly interacting electronic 
spin flip scattering (ESF) and the strongly interacting 
antiferromagnetic spin fluctuations (ASF) both of which 
show a frequency dependent Tj~^. The effects on both 
copper and oxygen [351 US] of ESF, where quasiparti- 
cles scatter onto the same node on the Fermi surface ef- 
fectively flipping spin but maintaining momentum, were 
found to have a uniform change of T^^ with tempera- 
ture. Fig. [sjb). Additionally, the Doppler shift is related 
to the angle between the scattering wave vector and the 
underlying crystal lattice giving a much larger distribu- 
tion in T^^. In the strongly interacting ASF limit, Fig. 
[Sf^a), the quasiparticles scatter between different nodes of 
the Fermi surface with a large change in momentum, and 
a sign change in the Doppler term of Eq. 1. This can give 
a non-zero value of relaxation at the minimum frequency 
point of the spectrum and a nonuniform change in T^^ 
with temperature. It was later determined |46| that, in 
the quantum limit, ESF can also cause a non-monotonic 
Tj~^ such that there is a local maximum in the rate at 



the minimum field point with the requirement that the 
line from the vortex along the node direction must cross 
the minimum of Hint- 



A different approach based on finding analytical and 
numerical solutions of the Bogoliubov-de Gennes equa- 
tions was taken by Throckmorton and VafekfTT^ calculat- 
ing the shifts and relaxation including strong magnetic 
fields up to 42 T. Their results show that in the near 
vortex core region antiferromagnetic fiuctuations are not 
necessary to account for experiment. i43) Rather they 
suggest that quasiparticle pair creation and annihilation 
(PGA) can be a dominant contribution to T^^ , especially 
important for high field NMR. 



In general the theory captures the experimental trends. 
However, there remain several questions at this point. Is 
the behavior of T^^ at the lowest frequencies in the spec- 
trum a manifestation of the ASF model? Gan the theory 
be extended to high magnetic fields and incorporate the 
Zeeman interaction for quasiparticles? And finally what 
is expected near the vortex core and is the PGA mecha- 
nism the dominant one? 
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FIG. 7. The magnetic field dependence of ^^O T^^ for two d- 
wave superconductors selectively measured at the saddlepoint 
in the magnetic field distribution where the Doppler terms 
mostly cancel by symmetry. Here there is evidence for a Zee- 
man contribution to the quasiparticle excitations which vary 
with magnetic field according to oc H^, reproduced from Oh et 
al. [44| The zero field limit of T^^ corresponds to thermal re- 
laxation and non-magnetic impurity effects, surprisingly simi- 
lar in these two different superconductors. The slopes for the 
two materials would be the same if the electronic g-factors 
were the same. 



Vortex Cores 

With the spatial resolution of T^^ established, the vor- 
tex core can be probed if the magnetic field is sufficiently 
large for NMR signal to be extracted from the nuclei that 
occupy the area within the vortex cores. The inter-vortex 
spacing, d, scales with magnetic field as o? = (f>Q/H. 
Taking the vortex core radius to be ^, the vortex core 
will occupy a fraction, ^ Tr^^/d^, of the cross-sectional 
area of the sample, which is the same fraction of the 
NMR spectrum that can be associated with nuclei in the 
core, the highest frequency part of the spectrum shown 
shaded in Fig. |6] This leads to a field dependent vor- 
tex fraction Ayortex / Atotai = 7ri?C^/(/)0j a quantity that 
grows linearly with external magnetic field. 

The spin 1/2 nucleus ^os^pj ^las been used by 
Kakuyanagai et al.^^ to probe the vortices in 
TlaBaaCuOe+s. The ^ostI nucleus, although not in the 
superconducting plane, is sensitive to antiferromagnetic 
fiuctuations from the ^"^Cu sites. By comparison, ^^O 
is less sensitive to antiferromagnetic fluctuations owing 
to cancellation of the singlet correlations between neigh- 
boring Cu. The 205^j NMR measurements of Tf^, per- 
formed in low magnetic field, H = 2.1 T, increased quite 
steeply with increasing frequency from the saddle point 
field. Spectra at different temperatures were discussed 
in terms of the full width at half maximum, Si, and the 
frequency width between points at 1% of the peak in 
the spectrum, Af. Below T = 20 K, Af becomes much 




48.8 48.85 48.90 48.95 49.00 49.05 49.10 49.15 



H|oc (in units of i^JX ) 




48.8 48.85 48.90 48.95 



.00 49.05 49.10 49.15 



H|oc (in units o^JX ) 

FIG. 8. The local magnetic field dependence of T^^ at dif- 
ferent temperatures T < Tm from Morr's theory. [45_ Bot- 
tom panel, intra-node scattering (electron spin flip, ESF). 
Top panel, inter-node scattering (antiferromagnetic spin-flip, 
ASF) processes. In the top panel, ASF scattering shows a 
non-monotonic increase of T^^ across the NMR spectrum 
causing a local minimum at the saddle point. In contrast the 
ESF scattering has a monotonic temperature dependence of 
Tj^^ such that it increases with temperature and frequency 
at all points in the spectrum. The curves are offset for clarity. 



broader than expected from a London model while Si 
changes rather little. Furthermore, T^^ is enhanced in 
the high frequency region with a peak at T = 20 K, taken 
to be the Neel temperature of the vortex core. 

Mitrovic et al.^\ have exploited availability of very 
high magnetic fields at the NHMFL to reach 42 T. These 
measurements have a substantial fraction of the ^^O 
NMR spectrum of YBa2Cu307_a; in the vortex core, 
nearly 16% of the total spectrum. The temperature de- 
pendence of {TiT)~^ at the saddle point were compared 
with that in the vortex core as shown at several magnetic 
fields in Fig. [9| Outside the core (TiT)~^ is a constant, 
but in the vortex core region, in magnetic fields above 

= 6 T, (TiT)-icx l/(r - e) following a Curie - Weiss 
law similar to ^'^Cu relaxation in the normal state, how- 
ever with < 0. This was interpreted as evidence of 
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FIG. 9. Tj^^ temperature dependence at different magnetic 
fields and vortex lattice locations. (TiT)"^ was found by 
Mitrovic et al. [33] to be a constant outside the vortex core in 
YBa2Cu307-;c. At = 6 T, (TiT)~^ is constant presum- 
ably because of insensitivity to the core at low magnetic field; 
however, with increasing H the temperature dependence fol- 
lows a 1/{T — 9), behavior of a Curie- Weiss law indicative of 
antiferromagnetism. 



strong antiferromagnetic fluctuations in the vortex core. 



vibrations show qualitatively the same behavior as Tj~ . 
However, the calculated rate is 3 orders of magnitude 
slower than what is experimentally observed. Further- 
more, the spin-echo decay rate follows neither a gaussian 
decay as expected for a ^^Cu T^^ mechanism [5Dj or an 
exponential decay due to vortex dynamics in the solid 
state. [2S] Although there is evidence that vortices play a 
role and that Tj"^ is spatially resolved, further work will 
be required to establish the basic mechanisms that are 
responsible for T2^- 

SDW in Bi2SrCa2Cu208+« 

The spatial dependence of ^^O relaxation in 
Bi2SrCa2Cu208+5 single crystals has been found to be 
significantly different than in YBa9 Cuq O7- ^ . [20] Both 
T]~^ and Tj"^ show a non-monotonic correspondence be- 
tween relaxation and local magnetic field, Fig. [lO) much 
more so than can be accounted for by any combination 
of Doppler shifts and broadened Redfield patterns for 
the distribution of local fields. It has been suggested 
that antiferromagnetic polarization near the vortex core 
could disturb the distribution of local magnetic fields and 
broaden the NMR spectrum. [T71 ST] Calculations along 
these lines were performed by these authors to account 
for broadening observed in the YBa2Cu307_2: ahgned 
powder experiments. According to one model. [T7] the 
magnetic moment in the vortex core required to account 



Spin-spin relaxation: spatial resolution 

Although the spin-spin relaxation, Tj"^, has not gained 
as much attention as T^^, it has also been shown to 
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have a spatial dependence. 



Specifically, T2 in 



both YBa2Cu307_^ and Bi2SrCa2Cu208+5 follow simi- 
lar frequency dependence at low temperat ure as T^^ . A 
direct comparison of Tj~^ and T2~^ in Fig. 10 emphasizes 
their similarity and naturally suggests that there are at 
least two significant contributions to Tjj"^, one of which 
has the same origin as for T^^ and which we believe 
depends on vortex supercurrents. The second is inde- 
pendent of frequency. The unusually pronounced upturn 
in the rates at low frequency in the NMR spectrum will 
be discussed in the next section. It was also shown that 
Tg"^ in YBa2Cu307_j; has similar behavior at different 
crystallographic positions, ^^Y and both ^""O sites, sug- 
gesting that the relaxation has a common origin likely 
determined by vortex vibrations. [39l |48] Although the 
theory[35] for vortex dynamics is not settled, as we shall 
discuss next, nonetheless its affect on Tj"^ is predicted to 
depend strongly on position in the vortex unit cell quite 
similar to calculations for T^^. 

Calculations fiU] based on Langevin dynamics of a fre- 
quency dependence of Tj"^ due to over damped vortex 




85.90 85.95 86.00 86.05 
Frequency (IVIHz) 

FIG. 10. The frequency dependent ^^O relaxation for 
Bi2SrCa2Cu2 08+a at T = 4 K and = 15 T from Mounce 
et aZ. [20] The spin-lattice relaxation rate (black, circles), T^^ 
and the spin-spin relaxation rate (blue, squares), T2^, have 
a strong non-monotonic relaxation as a function of fre- 
quency contrary to what is seen in YBa2Cu307_a; and 
Tl2Ba2Cu06+4. This behavior can be explained by a spin- 
density wave induced by the vortex core in its near neigh- 
borhood over a distance ~ 2^. In this model, the frequency 
distribution in the spectrum is relabeled non-monotonically, 
oscillating with position approaching the vortex core. 
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FIG. 11. The spectrum and frequency dependent T^^ relaxation rate at T = 4 K in Bi2SrCa2Cu2 08+i from Mounce et al. |20| 
for experiment (black) and their model calculation (red). With increasing external magnetic field, the ampHtude of the SDW 
component of the spectrum increases relative to the contribution from supercurrents. This leads to an increase of Tf^ at the 
lower frequency portion of the spectrum. 



for the broadening would need to be 2-5 /is if the anti- 
ferromagnetic wave vector is taken to be ^ tt/oo, where 
flo is the crystal lattice constant. However, the size of 
this moment is unreasonably large. [Sj 

There is independent evidence for vortices develop- 
ing spin-polarization around the vortex core. Indirectly, 
STM experiments in Bi2SrCa2Cu208-i-5 have shown a 
"checkerboard" pattern of the local density of states 
with periodicity 4ao which should be associated with an 
8ao spin density modulation. [511155] And there is direct 
evidence from neutron, elastic scattering experiments, 
performed on (Lai„a;Sra;)2Cu04 which indicate mag- 
netic Bragg scattering with a wavevector of 7r/8ao.[54]- 
I56j Furthermore, there have been theoretical predictions 
with models involving magnetic competition and/or co- 
existence with superconductivity consistent with these 
results. [57] 

The magnetic field dependence of the spectral 
linewidth at T = 4 K in clean single crystals of 
Bi2SrCa2Cu208+5 is well beyond what might be ex- 
pected from a GL-calculation of the NMR spectrum. A 
phenomenological model of the local magnetic fields was 
suggested by Mounce et o/.,[20j including contributions 
from vortex supercurrents superposed on a damped spin 
density modulation centered at the vortex cores, 

B{x, y) = Acos(27ra;/A)cos(27rt//A)e"(="'+2'')/'"' (3) 



Here A is the amplitude of the paramagnetic contribu- 
tion, A is the periodicity of the modulation and tr is 
the characteristic decay length away from the vortex 
core. Fitting these to the distribution of local magnetic 
fields results in A increasing approximately linearly with 
magnetic field while a is roughly constant and given by 
~ 2^. However, the addition of an oscillating local mag- 
netic field, Eq. 3, dramatically changes the dependence of 
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T-f^on frequency throughout the NMR spectrum, Fig 
and [TT] At low applied fields the relaxation profile is 
similar to that found in YBa2Cu307_2: as there is a cor- 
respondingly small contribution from vortex core, spin 
polarization. As external magnetic field increases, the 
spin density wave amplitude increases while the vortex 
supercurrent component decreases producing the non- 
monotonic behavior. In this work the relaxation rate was 
attributed to Doppler shifts of the nodal quasiparticles 
and was taken to be proportional to the square of the 
supercurrent momentum in the Cu02 plane, calculated 
as in Fig. [2] 



s-wave superconductors 

Without nodes in an energy gap at the Fermi surface, 
the quasiparticle bound states in an s-wave superconduc- 
tor are localized in the vortex core. The energy of the 
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metric with a constant T-^ throughout the spectrum, 



Fig. 12 



16.88 16.90 16.92 16.94 
Frequency (MHz) 



16.96 16.98 



FIG. 12. The NMR spectrum and relaxation, Tj"^ for the 
skutterudite compound, LaRu4Pi2 in relatively high magnetic 
field H = 0.98 T, H/Hc2 = 0.31, reported by Nakai et a/. ,60] 



core states were calculated to be E^^ = fiA^/Ep, with 
fi— 1/2, 3/2... , relative to the Fermi energy and are im- 
portant at low temperatures where they are evident in 
the zero-bias anomaly in tunneling experiments. |58L I59j 
The spatial dependence of Tj~^ would be homogeneous 
with the exception of the vortex core unless spin diffusion 
were to play a particularly important role; see Appendix 
B. 

Recent NMR in the s-wave superconducting skutteru- 
dite LaRu4Pi2 {Tc — 7.2 K) has a frequency dependent 
relaxation. [60] Above Tc the '^^P NMR spectrum is sym- 
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FIG. 13. The frequency dependent ^^As, relaxation for 
Bao.67Ko.33Fe2As2, at H = 16.5 T from Oh et al.^ At T = 4 
K, the spin-lattice relaxation rate has a strong frequency de- 
pendence across the NMR spectrum varying by one order of 
magnitude increasing with increasing frequency. This profile 
develops only in the superconducting state where = 38 K. 



As it is cooled below Tc the spectrum be- 
comes asymmetric with the familiar Redfield-like pattern, 
and T^^ develops a profile with a monotonic increase at 
higher frequencies. Since there are no nodal quasiparti- 
cles, the distribution of T^^ was attributed to diffusion 
from the vortex core. This is reasonable according to 
the theory of spin diffusion since these effects are larger 
for high gamma nuclei varying as the local dipolar field 
squared, ^ 7^. For ■^^P this is about an order of mag- 
nitude larger than it would be for ^^O. The experiments 
are done in relatively high magnetic field, H/Hc2 = 0.31 
where the vortex cores strongly overlap and half of the 
probe nuclei are within two coherence lengths of the core. 

The recently discovered pnictide high Tc 
superconductors [61] are thought to have s± super- 
conducting gap symmetry, a generically multiband 
superconductor with two principal gaps, one with hole 
excitations and the other around an electron pocket. [55] 
A calculation by Bang[53] indicates that the density 
of states will produce a Volovik-like effect where the 
quasiparticles have a finite density of states outside the 
vortex cores in a region where the supercurrents are 
sufficiently high to destabilize Cooper pairing from the 
smaller gap. It is especially important to include the 
interplay between interband scattering from impurities, 
temperature and magnetic fields. Bang's theory finds 
that the density of states has a linear field dependence 
at low magnetic fields, which a naive interpretation 
would suggest must produce a quadratic field depen- 
dence for the spatially averaged, rf^oc H^. However, 
spin-lattice relaxation in the mixed state in a multiband 
superconductor is not so straightforward and a more 
careful analysis[51] leads to a linear field dependence. 
Recent measurements on Bao.67Ko.33Fe2As2 of T^^ as 
a function of local magnetic field have revealed [H] an 
inhomogeneous distribution of relaxation that appears 
only in the superconducting state and this behavior is 
quadratic in field. In Fig. [13] we show the evolution of 
the Tf ^ profiles as a function of temperature providing 
evidence that Doppler effects play a role. However, 
discrepancy with the theory [64] is not understood. 



CONCLUSION 

This brief survey focuses on NMR in high tempera- 
ture superconductors emphasizing the effects of vortex 
supercurrents on the spectrum and their Doppler effects 
on quasiparticle excitations and correspondingly on 
spin-lattice and spin-spin relaxation. At this stage 
of development in the field it is clear that frequency 
resolved relaxation spectroscopy has important contri- 
butions to make to understanding vortex structures in 
unconventional superconductors. Frequency dependence 
implies spatial organization since the vortex currents 
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increase in amplitude approaching the vortex core oc 1/r, 
a consequence of the London equations. At the same 
time it is also clear that interpretation of experiments 
is greatly facilitated if samples are high quality single 
crystals with well-established electronic properties. 
Measurements at very high magnetic fields have been 
valuable and will continue to be significant particularly 
to improve our knowledge of vortex core excitations 
with a spin-dependent probe. In this regard NMR has 
a unique role to play. Our review discusses the subject 
of NMR investigations of vortices, mainly from the 
authors' perspectives and their contributions to it. 
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APPENDICES 
Volovik Effect 

For a simple metal T^^ is expressed in terms of the 
thermal average of the joint density of states: 

1/Ti (X J iV(e,)iV(6/)/(e,) [1 - f{ef)] de (4) 

where i and / indices label the initial and final electron 
states, /(e) is the Fcrmi-Dirac distribution function, and 
« e/. 



rate obtained from Eq. 4, i.e. < 1/Ti >av 

< l/Ti (X J < Nie^)N{e^) /(e^)[l-/(e^)] de 

(5) 

where the integral is a thermal average over all possi- 
ble quasiparticle states. Therefore we focus on just the 
important part of the integrand, 

< iV(e^)iV(eJ (X }^ j\N[e{r)]frdr (6) 

where e(r) is the spatial dependence of the energy in 
the vortex unit cell of size d and the integral goes from 
the vortex core radius, ^, to d. If the temperature is 
sufficiently low then the relevant quasiparticle excitations 
are small, we can expand A^(e) near e = 0: N{€) ~ e", 
provided that the density of states itself is not modified 
approaching a vortex core. For example, in the case of 
a clean d-wave superconductor the nodal quasiparticles 
have a dispersion, n — \. 

Now we introduce the inhomogeneity of the mixed 
state owing to the Doppler shift of the elementary ex- 
citations, 5e ^ • Ps, where oc 1/r. 

<N{e^)N{e^)>a. cx ^j\-^^rdr (7) 

For a d-wave superconductor n — 1, i.e. N{e) cx e, 
and the spatially averaged rate is given in Eq. 8 where 
we note that Hd^ = 4>q, and 4>q is the flux quantum: 

<l/Ti>av-d OC l^ln-^ oc H £n^ « H. (8) 

This is the accepted behavior for the spin-lattice relax- 
ation rate for a d-wave superconductor and was discussed 
by Volovik [33] for the spatially averaged density of states 
in the mixed state. 

Also for the d-wave superconductor we can write the 
average density of states as: 

<N[e)>a-o-d (X I, I r-^rdr ^ ]:^{d-0 (9) 



In a spin-lattice relaxation experiment we measure the 
NMR nuclear magnetization M{t) as a function of the 
time during relaxation after an i?_F'-excitation that dis- 
turbs M from equilibrium. Since the mixed state of a 
superconductor is spatially inhomogeneous we are con- 
cerned with the spatially averaged NMR response. The 
resulting M{t) is the superposition of the responses of 
nuclei from all regions of the sample. The initial time 
evolution of the magnetization, i.e. the limit as < — ^ 
of M(i), is equivalent to the magnetization evolution at 
the average rate. This allows us to compare experiment 
with a theoretical calculation of the spatially averaged 



<iV(e)>a„_d ~ J « H^''' (10) 

as expected. However, the approximations above 
and the conclusions reached are not valid for a si- 
superconductor. [M] 

On a historical note, a calculation along these lines was 
first made by Volovik[65 in 1988 in the context of en- 
couraging experimental efforts to explore the symmetry 
of the order parameter of the heavy fermion compound 
UPta . The Doppler effect from vortex supcrcurrents was 
then discussed theoretically by Yip and Sauls ,32] in 1992 
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for the d-wave case and then again by VolovikfS^ for the 
d-wave superconductor in his now famous paper in 1993. 

Finally, to emphasize a technical point, the average of 
T^^ is not equal to a rate determined from a standard 
relaxation experiment where the entire NMR spectrum is 
excited and detected. In fact, only the initial part of the 
recovery of the relaxation profile will correspond to the 
average rate, as we noted above. Alternatively, one can 
measure the frequency resolved rates across the spectrum 
and then calculate < T^^ >av as a spectrum- weighted 
average. 

Spin Diffusion 



lar field and AH is the maximal variation in the local 
field in the mixed state. Consequently, these authors 
ruled out the diffusion process in the earlier interpre- 
tations of the T^^ experiments. [661 [57] More recently 
WortisjTD] studied the possible effects of spin diffusion on 
T^^ in cuprate superconductors and compared her the- 
oretical calculations with experiment. [55 She considered 
YBa2Cu307_2; with realistic experimental conditions, al- 
lowing for spin diffusion coefficients determined from a 
combination of direct dipole-dipole coupling and indirect 
interactions. Her conclusions were the same as Genack 
and Redfield:[551 HHj «-e. spin difi'usion does not play a 
significant role in HTS spin-lattice relaxation. 



NMR measurements on superconductors and their in- 
terpretation include early work on the mixed state of low 
temperature superconductors, notably vanadium com- 
pounds, where T^^ was thought to be inhomogeneously 
distributed owing to spin diffusion. |661 167j The argument 
is that localized quasiparticle excitations in the vortex 
core, i.e bound states, [SB] would be a source of relax- 
ation which, through nuclear spin diffusion from outside 
the core region, would create a distribution in the nuclear 
Zeeman temperature. This effect would compromise the 
spatial identity of the relaxation profile discussed in this 
review of HTS. 

Diffusion can be attributed to the nuclear dipole-dipole 
interaction, via a 'fiip-flop' process, or similarly by any 
indirect electronic interaction, as is well known to ex- 
ist in some metals such as platinum and thallium. The 
measurements on vanadium superconducting alloys by 
Silbernagel et aL |66l I67| indicated an extra relaxation 
mechanism for T^^ at low temperatures beyond that ex- 
pected for an s-wave BCS state. This was qualitatively 
identified with nuclear spin diffusion of magnetization be- 
tween the vortex core and the bulk of the superconduc- 
tor. In these experiments a rough estimate was made of 
the diffusion coefficient from Tj"^ data. Later Genack 
and Redfield^Hl US] used field cycling techniques with 
samples of vanadium metal, exploiting the large gradi- 
ents that can exist in the mixed state of a type II su- 
perconductor. They measured the diffusion coefficient, 
D, and found reasonable agreement with the theory. But 
more importantly, they determined that spin diffusion in 
the superconducting mixed state is thermodynamically 
quenched in a very short time owing to depletion of the 
dipole energy reservoir. In these circumstances diffusion 
takes place with a diffusion constant, I?e// that is re- 
duced by a very substantial factor from its normal value, 

Deff - D (11) 

This reduction factor can be of order ^ 10^"^, and is pro- 
portional to the square of the product of natural abun- 
dance and the gyromagnetic ratio. Here Hd is the dipo- 
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